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Abstract 



In the companion paper, a coherent charge-current distribution for ra- 
diating electromagnetic wavelets was constructed on an oblate spher- 
oidal surface S a . Its main drawback was the necessity of including 
magnetic along with electric charges, making the sources impossible 
to realize. Here we show how this difficulty can be overcome by using 
Hertz potentials to generate a charge-current distribution due solely 
to bound electric charges. However, this distribution still appears dif- 
ficult to realize because it consists of multiple surface layers on S a . 
We show how it can be replaced by a simple volume distribution on 
a spheroidal shell. Our method generalizes the usual construction of 
equivalent Huygens sources, based on boundary conditions on an in- 
terface between electromagnetic media, by allowing the transition to 
be gradual without incurring addition complexity. 

1 Review of scalar wavelets 

We rely on the concepts in [K4], with some improvements in the notation. 
For further background on physical wavelets and complex-source pulsed 
beams, see [K94, HLKO, HF1]. The complex distance from the imaginary 
source point ia to the real observation point r will be denoted by 





and the complex time by 



The imaginary time b plays the role of an overall scaJe parameter, similar to 
the scale of ordinary wavelets in one dimension, determining the duration of 
the pulsed-beam wavelets. The imaginary space vector a similarly controls 
the spatial extent and orientation of the wavelets. The real and imaginary 
parts of r satisfy the inequalities 

\p\ < r, \q\ < a, where r = \r\, a = |a|. (3) 

For fixed the branch points of f in M 3 form the circle 

C = {r : r = 0} = {r : a ■ r = 0, r = a} 

and the 'standard' branch of r , defined by p > 0, has for its branch cut the 
disk 

V = {r : p = 0} = {r : a ■ r = 0, r < a}. 

r is real-analytic in R 3 except for a jump discontinuity due to a sign reversal 
upon crossing V. Every other branch satisfying the positivity condition 

a f^r>0 

can be obtained by continuously deforming V to a membrane B with the 
same boundary, 

8B = dV = C. 
The associated branch of f is defined by 

if, r<£V B 
r B = p B - iq B = i _ (4) 

y-r, reV B 

where V B is the compact volume swept out by deforming T> to B. It follows 
[K4a] that f B is real-analytic in M 3 except for a sign reversal across B. 

The scalar wavelet 1 with branch cut B is defined by 

_ 9(t-f B ) 



1 Strictly speaking, the term 'wavelet' should be reserved for certain choices of g, as 
explained in [K4]. The scalar wavelet of order n is obtained with g as the n-th derivative 
of the Cauchy kernel C(t) = 1/2-Kit. Also, note that we use units in which the propagation 
speed is c = 1 . 
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where g is the 'analytic signal' associated to a driving signal g(t) exciting 
the source by 

which is indeed analytic in the complement of the support of g: 

^ = J (d t - id b )g(t -ib)=0 Vt £ supp get. (6) 
at* 2 

The significance of the extension parameter b can be understood by noting 
that the real and imaginary parts of g are 

6 r g(f)df i r (f-t) g (f)df 

$5 is a smoothed version of g with 6 as the scale or resolution parameter, 
while % is a smoothed version of the Hilbert transform of g, again with 
b as the scale parameter. Thus, time variations of order less than \b\ are 
suppressed in g(t — ib). 

Due to its denominator, is singular on C, where f B = 0, and discontinuous 
in the interior of B, where r B reverses sign. To avoid any further singularities, 
we want to ensure that the numerator git — r B ) is analytic in all of M 3 , and 
for this it suffices to have its argument bounded away from the real axis by 
a positive distance. Since 

t-f B = t-p B -i(b-q B ) (8) 

and \q B \ = \q\ < a, a necessary and sufficient condition is 

a<|6|. (9) 

This states that the imaginary space-time four- vector (a, b) is time-like, 
belonging to the future cone of space-time if a < b and the past cone if 
b < —a. The condition (9) will be assumed from now on, making $? B real- 
analytic in M 3 — B. 

The source Sg of ^ B is defined by applying the wave operator: 

4vrS 8 = D* s where □ = d 2 t - A. (10) 
S B can be easily shown to vanish wherever ty B is twice differentiable, hence 

r iB => S 8 = 0. (11) 
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To characterize the source on B, we must apply □ in a distributional sense 
[GS64]. Just as differentiating the Heaviside function gives the delta func- 
tion, differentiating a discontinuous function like ^ B in a distributional sense 
gives a single layer on the surface of discontinuity, represented by a delta 
function of a variable normal to that surface. Since S B is obtained by dif- 
ferentiating ^ B twice, it will consist of a combination of single and double 
layers on B [K4, K4a]. Moreover, these layers diverge on the boundary 
C = dB since \T/ B is singular there. This singularity will be tamed below by 
combining wavelets with different branch cuts. 

The variables (p, q) defined by (1), together with the azimuthal angle <j) 
about the a-axis, determine an oblate spheroidal coordinate system where 
the level surfaces of p are the spheroids 

V ^-^2+- = 1 ' P^ ( 12 ) 

P Z + CT p Z 

and the level surfaces of q are the orthogonal hyperboloids TL q . All these 
quadrics are con focal, having the circle C as their common focal set. This is 
depicted in Figure 1. As p — * 0, S p shrinks to a double cover of the disk T>. 




Figure 1: (p, q, <fi) form an oblate spheroidal coordinate system. 



When the source of ^ B is computed, it will be singular on C due to the 
singularity of there. In the case of a real point source, this corresponds 
to the singularity of 5(r) at r = 0. Recall that the latter can be regularized 
by replacing the origin with a sphere of small radius r = a, whence 5{r) 
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is replaced by a uniform distribution on the spherical surface. The delta 
function can then be defined in terms of the limit a — > 0. The equivalent 
procedure now is to replace the sphere by the oblate spheroid S a , which is 
defined by p = Re r = a > 0. But we can go a step further and represent 
S a as a sum of two branch cuts, something that cannot be done for a real 
point source since the deformation of a point is still a point. Thus, consider 
the branch cut 

B + = S+ U A a (13) 
consisting of union of the the upper hemispheroid 

S+ = {reS a : z>0} 

and the apron 

A a = {r : r ■ a = 0, a 2 < r 2 < a 2 + a 2 } 

connecting <S e + to the branch circle C. (The apron must be included so that 
<9£> + = C as required.) 




Figure 2: The upper hemispheroidal branch cut B£ with its apron A a . 
Similarly, let 

B~ = S~ U A a (14) 
be the union of the lower hemispheroid 

S- = {r£S a : z<0} 
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with A a . For simplicity, denote the complex distance with branch cut 
by f± instead of f B ± and the corresponding wavelet by \l/ ± . Let be the 
interiors of the upper and lower hemispheroids and V a be the interior of S a . 
According to (4), 

reV a + => ¥+(f,t) = tf(-f,t), ¥_(f,t) = ¥(f,t) 
r G => *+(f,i) = *(f,i), *_(r,t) = *(-f,i) 
riV a => tf + (f,t) = tf_(f,t) = tf(f,t). 

Consider the average of * ± , 

*A(r,f) = -{¥_(f,t) + tf + (f,*)}. (15) 

Then by the above, 

*A{r,t) = < (16) 
[^ 2 {r,t), r f V a 

where the internal field and the external field ^2 are 

*i(f,t) = -{*(f,t) + *(-f,t)} , ¥ 2 (r,t) = *(f,t) (17) 

and will, for later purposes, be regarded as functions on all of 1R 3 . It follows 
directly from the definition (5) of g that 

^.f) = i^/_>«*'{f4--fnb} 

1 r (18) 



2vri (i - i') 2 - r 2 ' 

which depends only on f 2 and is therefore independent of the choice of 
branch cut. Furthermore, since (9) ensures that 

\t-t' ±r\ > |Im (t-t'±r)\ = \b ± q\ > \b\ - a > 0, 

we have 

|(t-t') 2 -r 2 | > (|6| - a) 2 >0 Vr. 

This shows that ^1 is real-analytic in R 3 , at least if g has compact support 
or decays sufficiently rapidly to ensure that the integral (18) converges. That 
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is, by taking the average (17) we have managed to cancel the singularities of 
*Sf(±r,t) on C as well as their jump discontinuities across V, leaving a field 
which is analytic in all of IR 3 and hence sourceless: 



□*i(f,t) =0V(r,t) eM 4 . 



The only region where ^a fails to be analytic is therefore S a , where it is 
discontinuous by (16). But even there, the irregularity is mild in the sense 
that the jump discontinuity 



is bounded — unlike that in any single branch *S> B , which is singular on C. 
Like and \I , 2, V^j will be regarded as a field on all of M 3 , although for the 
present we need it only on S a . The source T*a of ^a, defined by 



is therefore a distribution supported on S a . By the argument below (10), 
it consists of a combination of single and double layers on S a , with the 
difference that these layers are now bounded since S a avoids the singular 
circle C. Even so, it is not clear that the double layer can be realized in 
practice. In the next section we replace ^>a by a continuous field, where the 
transition from to ^2 occurs gradually over a range of spheroids, whose 
source is supported on a spheroidal shell instead of a single spheroid. (This 
explains the reason for viewing ^2 and * j as fields on all of ]R 3 .) Such 
volume sources, and their electromagnetic counterparts considered in the 
following section, should be realizable. 

2 Spheroidal shell sources 

Let H be the Heaviside step function. Since < p < a in the interior of S a 
and p > a in the exterior, we have 



It is natural to use the vector field orthogonal to S p given by [K3, Appendix] 



*j(r,t) = * 2 (r,i) - *i(f,t) = - {*(f,t) - *(-f,t)} 



(19) 



47rxu = , 



(20) 



^ A = H(a - p)^i + H (p - a)^ 2 - 



(21) 



n = Vp = 



pr + qa 



(22) 
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which is unnormalized with 

n 2 = \n\ 2 = P l + a l > 1, n-Vq = 0, Vn = ? 2p (23) 
jr + jr + <7 Z 

We can compute the source in (20) by using 

VH(p -a) = -VH{a - p) = H'(p - a)n = 5{p - a)n. (24) 

Applying the wave operator gives as a combination of terms involving 
6(p— a), interpreted as single layers on «S a , and 5'(p — a), interpreted as dou- 
ble layers. As mentioned above, it is doubtful whether the double layer can 
be realized in practice, and this will get still worse in the electromagnetic 
case, where the currents involve one more derivative. Since we are inter- 
ested in physically realizable sources, we now proceed to modify the above 
construction. The terms involving S(p — a) and S'(p — a) are unavoidable 
as long as the source is confined to the surface S a . To construct more real- 
istic sources, we now choose a function h that approximates the Heaviside 
function. A convenient example is 

h(p) = — &rg(— p + is) = —lm Log (—p + ie), e>0 (25) 

7T 7T 

which becomes the heaviside function in the limit as e — > 0. Note that 

is indeed an approximation to the delta function. But the simple choice 
(25) has the drawback that h! is not compactly supported, so the resulting 
source, although extremely small outside a shell of thickness O(e), will not 
have strictly compact spatial support. To obtain a compactly supported 
source, we now assume that h has the following properties: 



h(p) + h(-p) = 1 and h(p) = { ' 1 ~ (27) 




where e will be assumed fixed with 

< e < a, (28) 

so that h(p — a) vanishes on the disk V where p = 0. Define the regularized 
version of ^>a by replacing H by h in (21), 

m £ A = h(a-p)^i{r,t) + h{p-a)^ 2 {r r ,t). (29) 
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To simplify the equations, we use the abbreviations 

hi(p) = h(a -p), h 2 (p) = h{p - a), (30) 

so that 

V £ A = h^! + h 2 ^ 2 = h k V k (31) 

where the Einstein convention of summing over repeated indices has been 
used. The source T> A of ^> e A is defined as usual by 

4ttE^ = nV e A . (32) 

To compute this, note that since (27) implies h'(—p) = h'(p), 

= -h'(a - p)^in + h'ip - a)^ 2 n + h k V^ k 

= h'(p - a)^jn + h k VV k (33) 

where defined as in (19) and given by 

*j(f, t) =-L. g{t') dt' { ^ + 1 

47rw./-oo {t-t'-f t-t' + fj 

_ J_ [°° {t-t')g{t')dt' 

~ 2 7 rzf7_ 00 {i-t'f-r 2 ' 1 ' 

no longer represents a jump discontinuity of the field since we are not con- 
fined to a single spheroid. By the same argument used to show that is 
real-analytic in IR 3 , it follows that ^ j is real-analytic in IR 3 except for being 
discontinuous on V and singular on C due to the factor 1/f. Taking the 
divergence of (33) gives 

A^ £ A = h"ip - a)^jn 2 + 2h'(p - a)V*j • n + h'ip - a)fjV • n + h k A^ k . 
Since ^j(f, t) is complex-analytic in f when r ^ P, 

p > => V*j = *j Vf 
where the prime denotes the complex derivative with respect to f, 

*/ = ^ = ^ + id q )9j = \ {*'(f, f) + tf(-f , *)} . (35) 

By (23), 

Vf • n = (Vp — iVg) • n = n 2 . 
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Subtracting df^ e A thus gives 

-4tt£^ = h"(p - a)^jn 2 + 2h'(p - a)^/n 2 + h'(p - a)* jV ■ n - h k U^ k . 

But we have seen that vanishes identically and is supported on 
V, where h 2 = by (27) and (28). Using (23) gives the regularized source 



-4tt£ £ a = h"(p - a)*jn 2 + 2h ] P t 2 Q) {(p 2 + a 2 )*/ + ptfj} 

p z + q z 



(36) 



supported on the spheroidal shell 

5^ = {r : q - e < p < a + e}. (37) 

We emphasize that T, A is a smooth volume source that depends only on the 
'jump field' tyj. Taking the limit e — > so that /i becomes the Heaviside 
function gives the source T, A consisting of single and double layers on S a . 

3 Maxwell's equations and Hertz potentials 

We work with the following complex combinations of electromagnetic fields: 

F = D + iB (38) 
G = E + iH = F- AttP (39) 
P = P e + iP m (40) 

where the units are Gaussian with c = 1, P e is the electric dipole density, 
and P m is the magnetic dipole density. Maxwell's equations take the form 

V • F = 4irp F + iV xG = -4vr J (41) 

where F = dtF. In the general case of complex charge and current densities 

p = p e + ip m J = J e + iJ m , 

equations (41) are equivalent to 

V • D = 4irp e D-V x H = -4irJ e (42) 

V • B = Airp m B + V xE = -4vr J m , 
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so the imaginary parts (p m , J m ) represent the magnetic charge-current den- 
sity. Since magnetic monopoles are not observed, we must require 

p m = and J m = 0. (43) 

That is, (41) are completely equivalent to the usual Maxwell equations if we 
add the requirement that (p, J) is real. We will consider solutions derived 
from a complex Hertz potential consisting of electric and magnetic Hertz 
vectors 

Z = Z e + iZ m (44) 

whose source is the polarization, 

DZ = 4vrP. (45) 

(For this reason, (Z e ,Z m ) are sometimes called polarization potentials.) 
The field F is then given in terms of Z by 

F = VxVxZ-HVxZ, (46) 

and it follows from (39) and (45) that 

G = F-DZ = VxVxZ + iVxZ + AZ-Z 

= VV • Z + iV x Z - Z. (47) 

The real form of equations (45), (46) and (47) is [BW99, pp 84-85] 

DZ e = 4vrP e UZ m = 47rP m 

E = VV-Z e -VxZ m -Z e B = VxVxZ m + VxZ e 

H = VV • Z e + V x Z e - Z m D = V xV x Z e -V x Z m . 

An inspection of the expressions for E and B reveals the meaning of the 
Hertz potentials as 'superpotentials' from which the four-vector potential 
($>, A) can be derived by 

$ = -V • Z e A = V x Z m + Z e . (48) 

In fact, these automatically satisfy the Lorenz condition 

6 + V- A = 0, 
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and every four-vector potential satisfying it can be derived from Hertz po- 
tentials. The freedom to choose a gauge for (<E>, A), including a non-Lorenz 
gauge, is part of a much greater gauge freedom in (Z e , Z m ) [N55, BW99]. 

According to (46), F is a curl, so by (41) the free charge density vanishes: 

p = 0. (49) 

Furthermore, (47) gives 

VxG = iVxVxZ-VxZ = iF, 

therefore by (41), the free current density also vanishes: 

J = 0. (50) 

Maxwell's equations (42), written in terms of the microscopic fields (E, B), 
now state that 



V -B = 

V ■ E = 4irp b 

where 

Pb = -V ■ P e 



V x E + B = 

V xB-E = 4tt J b (51) 
J b = P e + V x P m (52) 



represent the bound charge and current densities generated by the variable 
polarizations (P e ,P m ). The fields derived from Hertz potentials as above 
are thus due entirely to bound sources. 2 



4 Spheroidal electromagnetic antennas 

In this section we construct electromagnetic wavelets from scalar wavelets by 
turning ^ into Z, then compute their charge-current densities. It is essential 
that the polarization P defined in (45) have compact spatial support, as it 
can otherwise not be realized. There are various ways to turn a scalar 
solution of the wave equation into a vector solution without increasing the 
support of its source distribution, the simplest being 

Z = p$> (53) 
2 Free charge-current densities can be added by using stream potentials [N55]. 
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where p is a constant (possibly complex) vector. The polarization is then 
given by 

4vrP = pU$ = 47rpE, (54) 

so P and £ have the same support. Since S is a distribution consisting 
of single and double layers on T>, so is P. A similar construction applies 
to the different versions supported on the general branch cut B and the 
spheroid <S Q . As explained below (24), the layers on B are singular on C 
while those on S a are bounded. Even so, the charge-current distributions 
(52) require one further differentiation, hence they generate a still higher 
layer with coefficient S"(p — a), and it is doubtful whether such distributions 
can be realized. For this reason we confine our analysis to volume sources 
on the spheroidal shell (37). Define the Hertz potential 

Z £ A = P n (55) 
with ^ A as in (29), whose polarization density is 

P £ a=P^a (56) 

with T, A given by (36). If we interpret T, A as a scalar density, then (56) 
suggests an interpretation of p as a (complex) combination of electric and 
magnetic dipole moments. The charge and current densities on the shell, as 
given by (52), are 

pb = _ Re {V • P A } = - Re {p ■ VS^} (57) 

and 

J b = Re{p^} + Im{V x P £ A } = RejpS^} - Im{p x VS^}. (58) 

Outside the shell the potential Z £ A coincides with Z = p^, whose pulsed- 
beam field F was computed in [K4]. 

5 Extended Huygens sources 

The above suggests an generalization of Huygens sources [HY99], allowing 
equivalent sources to be represented on shells instead of surfaces surround- 
ing a bounded source. We present this generalization and compare it to the 
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usual method based on boundary conditions on an interface between electro- 
magnetic media. Let p(r, t) be a differentiable function, which will be called 
a zone function. Fix two numbers p\ < P2 and consider the time-dependent 
surfaces and volumes in M 3 defined by 

5i(t) = {r : p(r, t) = Pl }, S 2 (t) = {r : p(r, t) = p 2 } 
Vi(t) = {r : p(r, t) < Pl }, V 2 (t) = {r : p(r, t) > p 2 }. 

Given two electromagnetic fields (Fi,G\) and (F 2 ,G 2 ), with or without 
sources, we want to construct an interpolated field (F, G) so that 



F(r,t) 



Fi(r,i), re^i(t) 
F 2 (r,t), reV 2 (t) 



G(r,t) = \ G ^ T ^ (59) 
\G 2 {r,t), reV 2 (t). 



Choose a differentiable function h 2 (r,t) such that 



"™ (60) 
[1, r£V 2 (t) 

and let 

/ii(r,t) = l-h 2 (r,t). 
We define the interpolated field as 

F(r,t) = h k (r,t)F k (r,t) 

G{r,t) = h k (r,t)G k (r,t) (61) 
where summations over k = 1, 2 are implied, and the jump field 

Fj = F 2 -F 1 = D J + iBj 

Gj = G 2 -G 1 = E J + iHj. (62) 
Then, according to (41), the charge density of (F,G) is 

Airp = V • F = 4irh kPk + Vh 2 ■ Fj (63) 

where 

A7r Pk = V-F k , k = 1,2 (64) 
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are the charge densities of the prescribed fields. Thus, in addition to the 
interpolated charge density 

pi = h k pk (65) 

we have a transitional charge density given by 

Airpr = Vh 2 ■ Fj (66) 

which depends only on the component of the jump field F j parallel to Vh 2 . 
According to (60), px vanishes outside the transition shell 

V T (t) = {r: Pl <p(r,t)<p 2 }. (67) 

Similarly, the current density is 

4ttJ = -F - iV x G = 4irh k J k - h 2 F j - iVh 2 x Gj (68) 

where 

4vrJ fc = -F k - iV x G k , fc = l,2 

are the current densities of the prescribed fields. Hence J is the sum of the 
interpolated current density 

Ji = h k J k (69) 

and a transitional current density on Vp (t) given by 

4tt J t = -h 2 Fj - iVh 2 x Gj (70) 

which depends only on Fj (if h 2 is time-dependent) and the component of 
Gj orthogonal to Vh 2 . The electric and magnetic transitional sources are 
obtained by taking real and imaginary parts. Assuming h 2 is reaJ, this gives 

47T/4 = Vh 2 ■ Dj inJ? = -h 2 Dj + V/i 2 x Hj 

inp^p = Vh 2 ■ Bj 4irJ% = -h 2 Bj - Vh 2 x Ej. (71) 

Letting h 2 be complex in Vr(t) makes the transition shell a chiral medium 
mixing electric and magnetic fields. A further generalization is obtained by 
replacing h k with 3x3 matrices (dyadics) M k satisfying 

U 2 (r,t) = r 'I*®, H 1 (r,t)=I-H 2 (r,t) (72) 
II, re V 2 (t) 
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where I is the unit matrix. This makes the transition shell Vr{t) a non- 
isotropic medium as well as chiral if are complex. See [LSTV94] for a 
treatment of chiral and non-isotropic media. 

Choosing the zone function p(r,t) time-dependent thus gives a simple for- 
mulation of the transition shell as a moving source, which could be useful 
in the analysis of radiation by moving objects. 

To see how all this relates to Huy gens' principle, suppose we are only given a 
field (F2, G2) whose charge-current density (p2, J 2) is confined to V±(t), and 
want to find an equivalent charge-current density confined to Vr(t) whose 
radiated field in V^i) (but not necessarily elsewhere) is (F21G2). We are 
free to choose the field (Fi,Gi) in any way that gives vanishing interpolated 
sources 



since the sources of the interpolated field are then purely transitional and 
hence confined to Vr(t) as desired. To satisfy (73), it suffices to require that 
(pi, Ji) be confined to V^i). Thus, choosing any field (Fi,G±) with sources 
in V2(t) and any function h,2 satisfying (60), an equivalent charge-current 
density on Vr(t) is given by (66) and (70). As p\ — ► P2, Vr(t) becomes S2 
and (pt,Jt) become ordinary Huygens surface sources [HY99]. 

The freedom to choose (F±,Gi) (interpreted as the 'interior field' if Vi(t) 
is bounded) and I12 is constrained by the requirement that the magnetic 
charge-current density must vanish, as detailed below. 

Now suppose that p = p(r) is time-independent, so Sk and Vk are fixed, and 
choose hk to be time- independent and real. As p\ — > P2, assume that 



where n(r) is a normal vector field on 52 pointing into V2. Then (66) and 



Pi = h k pk = 0, 



J 1 = hkJk = 



(73) 



lim Vh2(r) = 5(p(r) — P2)n(r) 



(70) give 



PT -> S(p(r) -p 2 )cr 
J T -> S(p(r) -p 2 )K, 



(74) 



where 



Alio = n ■ Fj and 4ttK = —in x Gj. 



(75) 
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are the surface charge and current densities on S 2 , whose real and imaginary- 
parts give the electric and magnetic surface sources: 

Aira e = n ■ Dj 47rK e = n x Hj 

Aira m = n B,j 4vrK m = -n x Ej. (76) 

Since magnetic monopoles are not observed, a m and K m must vanish. When 
Pi < P2 , this may be accomplished if h 2 can be chosen so that 

Vh 2 -Bj = 0, Vh 2 xEj = 0, (77) 

which is possible 3 if 

Ej-Bj = Vr£V T (t). (78) 

In the limit p\ — > p 2 , (77) reduces (76) to 

n ■ Dj = 4Tra e n x Hj = 4irK e 

n-Bj = nx Ej = 0, (79) 

which are the usual boundary conditions on an interface between two media. 

Returning to the general time-dependent setting, consider now an alternative 
procedure of special interest here. Instead of interpolating two prescribed 
fields (Fk,G}~), let us interpolate two Hertz potentials Z^: 

Z = h k Z k . (80) 

As seen, this automatically results in vanishing 'free' sources p = 0, J = 0. 
The polarization is found to be 

4vrP = DZ = 4Trh k P k + 2h 2 Zj - 2(V/i 2 ■ V)Zj - Z jAh 2 

where ^irP k = OZ k are the polarizations of the prescribed fields. If 

suppPi C V 2 and suppP 2 C V\, 

then the interpolated polarization h k P k vanishes and the polarization is 
purely transitional on Vr(t): 

4irP = 2h 2 Zj - 2(Vh 2 ■ V)Zj - ZjAh 2 . (81) 

This generalizes (56), and the bound charge-current densities derived from 
P e and P m via (52) generalize (57) and (58). 

Letting Vr be a nonisotropic medium by using Mk (72) makes it easier to enforce the 
absence of magnetic monopoles. 
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6 Conclusions 



We have improved on the computation of sources for electromagnetic wavelets 
given in [K4] in two ways: (a) The spheroidal surface S a supporting the 
sources has been replaced by a spheroidal shell <S„ supporting smooth vol- 
ume sources. This eliminates the multiple layers on S a which make the 
sources difficult if not impossible to realize, (b) By deriving the sources 
from Hertz potentials, we have eliminated the magnetic charge-current den- 
sity, further facilitating their realizability. 

The problem with the magnetic sources in [K4] can be better understood 
from the current perspective. Let the zone function be p = Re f and S k , V k 
be as above with < p\ < P2- Let the fields (F k , G k ) be derived from the 
Hertz potentials 

Z fc =ptf fc , k = l,2 

with ^ k given by (17). Recall that ^ k are analytic in (r,t) for p > 0. 
Therefore the fields (F k ,G k ) are analytic in (f,t) for r V, and 

P k = 0, G k = F k = E k + iB k Vr(£V. 

The jump fields 

Fj = Gj = Ej + iBj 

are also analytic, as is their polarization scalar [K3a] 

F 2 =Fj-Fj = E 2 -B 2 + 2iEj ■ Bj. (82) 

The condition (78) thus requires the imaginary part of an analytic function 
to vanish for r £ Vt, which implies that Fj vanishes identically. Although 
every electromagnetic field must have F 2 — > in the far zone [B15], fields 
satisfying F 2 = globally, called null fields, are rather degenerate. In 
particular, the electromagnetic wavelet fields and not null and hence cannot 
fulfill (78). Instead, we have begun with an interpolated Hertz potential (80) 
and derived bound sources in the transition shell, thus preserving analyticity 
without invoking the existence of magnetic monopoles. 

There is still an unsatisfactory aspect to the polarization (56) and charge- 
current density (57), (58). Namely, they depend on the fixed vector p and 
thus do not conform to the spheroidal geometry. This suggests using meth- 
ods of constructing Z from ^ other than (53). While (53) is the complex 
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version of Whittaker's potentials [W4], there are alternatives which do not 
require a fixed polarization vector, such as Debye potentials; see [BD98] for 
example. Such alternatives will be considered in future work. 

Note: After this paper was finished, I learned from Dr. Arthur Yaghjian that 
a similar formulation of the gradual transition between two electromagnetic 
media has been developed by Lindell, Tretyakov and Nikoskinen [LTNO]. 
The present form is somewhat more general in the following respects, (a) 
The regions and surfaces are allowed to be time-dependent, (b) The deriva- 
tion of generalized Huygens sources in [LTNO] makes the assumption that 
one of the prescribed fields vanishes, which we find to be unnecessary. This 
allows extended Huygens sources on Up with nonvanishing 'internal' fields. 
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